Numerical Methods for Multilattices

Among the efficient numerical methods based on atomistic models, the quasicontinuum (QC) method has attracted growing interest in recent years. The QC method was first developed for crystalline materials with Bravais lattice and was later extended to multilattices (Tadmor et al, 1999). Another existing numerical approach to modeling multilattices is homogenization. In the present paper we review the existing numerical methods for multilattices and propose another concurrent macro-to-micro method in the numerical homogenization framework. We give a unified mathematical formulation of the new and the existing methods and show their equivalence. We then consider extensions of the proposed method to time-dependent problems and to random materials.Comment: 31 page

Matching Conditions in Atomistic-Continuum Modeling of Materials

A new class of matching condition between the atomistic and continuum regions is presented for the multi-scale modeling of crystals. They ensure the accurate passage of large scale information between the atomistic and continuum regions and at the same time minimize the reflection of phonons at the interface. These matching conditions can be made adaptive if we choose appropriate weight functions. Applications to dislocation dynamics and friction between two-dimensional atomically flat crystal surfaces are described.Comment: 6 pages, 4 figure

Absorbing boundary conditions for the Westervelt equation

The focus of this work is on the construction of a family of nonlinear absorbing boundary conditions for the Westervelt equation in one and two space dimensions. The principal ingredient used in the design of such conditions is pseudo-differential calculus. This approach enables to develop high order boundary conditions in a consistent way which are typically more accurate than their low order analogs. Under the hypothesis of small initial data, we establish local well-posedness for the Westervelt equation with the absorbing boundary conditions. The performed numerical experiments illustrate the efficiency of the proposed boundary conditions for different regimes of wave propagation

Impact of Numerical Methods in Thermal Modeling of Li-Ion Batteries on Temperature Distribution and Computation Time

Thermal battery modeling is important for further battery development and optimization. The temperature strongly influences the performance and aging behavior. In the cell stack, electrochemical processes take place resulting in a large amount of heat release, which, in turn, affects the temperature distribution. Therefore, the main focus is on the cell stack, the most complex structure inside the cell. In particular, the discontinuous and anisotropic material properties represent a major challenge for simulations due to the layering. This work proposes self-developed methods, based on the Finite Volume Method and the Finite Element Method, taking on these challenges. First, for both methods the functionality is verified and numerical convergence is validated. These, and also classical methods, are compared based on test problems with a known analytical solution in view of numerical errors as well as computing time. It if found that their accuracy and efficiency depends strongly on the specific problem, which makes their numerical investigation necessary and inevitable. Second, the methods are evaluated on a specific battery problem. Their results are plausible and correspond to the physical phenomena

Propagating modes of non-Abelian tensor gauge field of second rank

In the recently proposed extension of the YM theory, non-Abelian tensor gauge field of the second rank is represented by a general tensor whose symmetric part describes the propagation of charged gauge boson of helicity two and its antisymmetric part - the helicity zero charged gauge boson. On the non-interacting level these polarizations are similar to the polarizations of the graviton and of the Abelian antisymmetric B field, but the interaction of these gauge bosons carrying non-commutative internal charges cannot be directly identified with the interaction of gravitons or B field. Our intention here is to illustrate this result from different perspectives which would include Bianchi identity for the corresponding field strength tensor and the analysis of the second-order partial differential equation which describes in this theory the propagation of non-Abelian tensor gauge field of the second rank.Comment: 22 pages, Latex fil

Projected SO(5) Hamiltonian for Cuprates and Its Applications

The projected SO(5) (pSO(5)) Hamiltonian incorporates the quantum spin and superconducting fluctuations of underdoped cuprates in terms of four bosons moving on a coarse grained lattice. A simple mean field approximation can explain some key feautures of the experimental phase diagram: (i) The Mott transition between antiferromagnet and superconductor, (ii) The increase of T_c and superfluid stiffness with hole concentration x and (iii) The increase of antiferromagnetic resonance energy as sqrt{x-x_c} in the superconducting phase. We apply this theory to explain the ``two gaps'' problem found in underdoped cuprate Superconductor-Normal- Superconductor junctions. In particular we explain the sharp subgap Andreev peaks of the differential resistance, as signatures of the antiferromagnetic resonance (the magnon mass gap). A critical test of this theory is proposed. The tunneling charge, as measured by shot noise, should change by increments of Delta Q= 2e at the Andreev peaks, rather than by Delta Q=e as in conventional superconductors.Comment: 3 EPS figure

Large spin limits of AdS/CFT and generalized Landau-Lifshitz equations

We consider AdS_5 x S^5 string states with several large angular momenta along AdS_5 and S^5 directions which are dual to single-trace Super-Yang-Mills (SYM) operators built out of chiral combinations of scalars and covariant derivatives. In particular, we focus on the SU(3) sector (with three spins in S^5) and the SL(2) sector (with one spin in AdS_5 and one in S^5), generalizing recent work hep-th/0311203 and hep-th/0403120 on the SU(2) sector with two spins in S^5. We show that, in the large spin limit and at the leading order in the effective coupling expansion, the string sigma model equations of motion reduce to matrix Landau-Lifshitz equations. We then demonstrate that the coherent-state expectation value of the one-loop SYM dilatation operator restricted to the corresponding sector of single trace operators is also effectively described by the same equations. This implies a universal leading order equivalence between string energies and SYM anomalous dimensions, as well as a matching of integrable structures. We also discuss the more general 5-spin sector and comment on SO(6) states dual to non-chiral scalar operators

Complexity Analysis of a Fast Directional Matrix-Vector Multiplication

We consider a fast, data-sparse directional method to realize matrix-vector products related to point evaluations of the Helmholtz kernel. The method is based on a hierarchical partitioning of the point sets and the matrix. The considered directional multi-level approximation of the Helmholtz kernel can be applied even on high-frequency levels efficiently. We provide a detailed analysis of the almost linear asymptotic complexity of the presented method. Our numerical experiments are in good agreement with the provided theory.Comment: 20 pages, 2 figures, 1 tabl

Circular and Folded Multi-Spin Strings in Spin Chain Sigma Models

From the SU(2) spin chain sigma model at the one-loop and two-loop orders we recover the classical circular string solution with two S^5 spins (J_1, J_2) in the AdS_5 x S^5 string theory. In the SL(2) sector of the one-loop spin chain sigma model we explicitly construct a solution which corresponds to the folded string solution with one AdS_5 spin S and one S^5 spin J. In the one-loop general sigma model we demonstrate that there exists a solution which reproduces the energy of the circular constant-radii string solution with three spins (S_1, S_2, J).Comment: 16 pages, LaTeX, no figure